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ABSTRACT 

The shape of a cusped galactic nucleus is constrained by the range of shapes of orbits in 
its gravitational potential. It is shown for scale-free non-axisymmetric discs that while 
a plausible elongated density model requires at least some orbits to spend more time 
near the major axis than anywhere else, both regular boxlets and tube orbits generally 
cross the major axis too fast for self-consistency. If galaxies host inner nuclear discs or 
flat bars with a cuspy surface light profile (7 = dlog/i/dlogr > 0), their ellipticity 
1 — b/a cannot be greater than about 7/2. Discs or bars with a shallow central profile 
(7 < 0.3) should not be strongly elliptical. 

Key words: galaxies: kinematics and dynamics 



1 INTRODUCTION 

Current formation theories emphasize the roles of dissipa- 
tion and galaxy interaction as major processes in shaping 
present day galaxies, but is there any stringent limit on the 
shapes and density profiles of galaxies from general condi- 
tions such as equilibrium and stability alone? Particularly, 
are there stable triaxial equilibria with realistic radial den- 
sity profiles? This has been an out-standing stellar dynami- 
cal question ever since Binney (1978) invoked triaxial equi- 
libria to account for the observed flattened shape of ellipti- 
cal galaxies and their lack of rotation. Theoretically triaxial 
equilibria exist at least for systems with a finite density core 
and a wealth of box orbits and tube orbits. This was demon- 
strated for general models by Schwarzschild (1979, 1982), 
and for models with separable potentials by Statler (1987). 
However, the traditional assumption of a finite core in ev- 
ery galaxy was challenged by recent observations of nuclei of 
nearby elliptical galaxies with the Hubble Space Telescope. 
It was found that giant ellipticals have a power-law surface 
density distribution [i oc r~ 7 with < 7 < 0.3 near the 
center; the slope steepens to as much as 0.3-1.3 for small 
ellipticals, with the dividing line at Mb ~ —20 mag (Crane 
et al. 1993; Jaffe et al. 1994; Lauer et al. 1995; Carollo et 
al. 1997; Faber et al. 1997). These observations call for a 
re-examination of the existence of triaxial equilibria in a po- 
tential with a divergent force (F oc /i — *■ 00 for 7 > 0) at 
the center (Gerhard & Binney 1985; de Zeeuw & Carollo 
1996). Whether triaxial models of this kind exist also be- 
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comes a key uncertainty in, e.g., interpreting the kinematic 
and photometric data of galactic nuclei, weighing their cen- 
tral black holes and reconstructing the formation history of 
these systems. 

The conclusion that at least some strongly non-axi- 
symmetric models with a steep central cusp are proba- 
bly not in rigorous steady state is based on a handful of 
three-dimensional dynamical models. These have been built 
with various implementations of Schwarzschild's method, 
in which individual orbits are populated so as to match 
the model density, and include three-dimensional scale-free 
models with logarithmic potentials (Schwarzschild 1993), 
and two non-scale- free models (Merritt & Fridman 1996). 
Unfortunately, the power of these few numerical experi- 
ments is limited when it comes to exploring the parameter 
space. It is not clear how to extrapolate results obtained 
for a few models to a general statement about the whole 
class of cusped triaxial potentials, because the meaning of a 
small mismatch in the reconstructed density, which is often 
of the order of one percent or less, has to be interpreted on a 
case-by-case basis. Kuijken (1993, hereafter K93) showed in 
his systematic study of two-dimensional non-axisymmetric 
models with a logarithmic potential that whether the nu- 
merically constructed model is in equilibrium depends sen- 
sitively on numerical details, including the resolution of the 
spatial grid for the mass model, the grid for orbital initial 
conditions, the number of orbits used, and the integration 
time for each orbit. Although it may be feasible with present- 
day computer technology to carry out a massive numerical 
search in the multi-dimensional parameter space (axis ratios, 
inner and outer density slopes with the possible addition of 
a central black hole and the tumbling speed of the poten- 
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tial), it is intrinsically difficult in this approach to pinpoint 
the exact origin of any mismatch between the orbits and the 
density model. For example it has not been well-understood 
why replacing box orbits in a cored potential by boxlets in 
a cusped potential upsets the equilibria. 

In this paper we present a new approach to study self- 
consistency of a general non-axisymmetric model. We re- 
strict ourselves to the two-dimensional case of scale-free 
discs, and show that a non-trivial and necessary condition 
for the existence of a self-consistent elongated disc is that 
the angular speed of the regular boxlet and tube orbits when 
they cross the major axis should be consistent with the lo- 
cal curvature of the density distribution. It is well known 
that boxlets are less useful than box orbits when it comes 
to fit the model density near the major axis because the 
boxlets have their density maxima at the turning points 
rather than on the major axis (e.g., Pfenniger & de Zeeuw 
1989; Schwarzschild 1993). K93 made the interesting obser- 
vation that a box orbit has only one 'corner' per quadrant, 
while a boxlet orbit has two or more correlated 'corners' per 
quadrant, which makes them less flexible in fitting the model 
density. He also suggested that the lack of self-consistency in 
elliptical models is likely due to the spiky angular distribu- 
tion of boxlets rather than to the lack of flattened boxlets. 
Syer & Zhao (1998, hereafter SZ98) suggested that the self- 
consistency is first broken down near the symmetry axes of 
the model. Unfortunately, the result of SZ98 is limited to a 
special subset of the separable non-axisymmetric potentials 
introduced by Sridhar & Touma (1997, hereafter ST97). 

The structure of this paper is as follows. §2 gives a rig- 
orous formulation of the problem for scale-free discs. §3 il- 
lustrates the requirement on the curvature of the model with 
specific orbits in elliptic disc potentials. §4 gives the results 
of fitting the curvature, and §5 examines the assumptions in 
the model, and discusses generalizations of the method to 
three-dimensional and non-scale-free systems, as well as the 
implications for barred galaxies. 



2 SCALE-FREE NON-AXISYMMETRIC DISCS 

2.1 Formulation 

Consider a general two-dimensional scale-free non-axisym- 
metric disc potential 



(r,0)ocr a p(0)> 



-1< a< 1, 



(1) 



where p(0) defines the angular shape, with = and 7r/2 
the directions of the minor and major axis, respectively. We 
consider surface densities of the form 



^(r,e)xr- J s(e), 



< 7 < 2, 



(2) 



where 7 is the cusp strength and s(0) the angular shape of 
the surface density. In self-consistent systems s(0) and p(9) 
are related, and 7 = 1 — a. 

The regular orbits in this potential can be grouped ac- 
cording to their shapes (bananas, fishes, tubes, etc.; see 
Miralda-Escude & Schwarzschild 1989; K93), irrespective 
of their sizes, into self-similar families. Each family can be 
characterized by a dimensionless second integral, say, I, and 
the whole family can be then built by rescaling one reference 
orbit with a trajectory described by the polar coordinates 



(77 (t), 9i(t)). The weights among the self-similar 'cousin' or- 
bits should be prescribed in a scale-free fashion such that 
each family produces a r a_1 power-law density distribution, 
as for the model density. The different families should be 
weighed according to a to-be-determined positive function 
W (I) > so as to reproduce the angular part of the model 
density. 

A clear discussion of the algorithm for constructing self- 
similar non-axisymmetric models was given by Richstone 
(1980) and Schwarzschild (1993) for three-dimensional mod- 
els, and by K93 for two-dimensional models. These previous 
studies were focused exclusively on the logarithmic poten- 
tial, and it appears that no one has applied the algorithm 
to other non-axisymmetric power-law potentials, which are 
most relevant for galactic nuclei. Following K93, we first di- 
vide the angular coordinate into n slices with an interval 
A9 — 2n/n, and n approaching infinity. We then compare 
the amount of mass a reference orbit deposits in each an- 
gular sector to the amount of mass required by the density 
model in the same sector within the radius of the orbit. If a 
reference orbit (ri(t),9i(t)) spends a fraction A(t/T)i.j of 
the integration time T in the j-th sector, then it deposits 
an amount proportional to w(I)A(t/T)i.j. By comparison, 
the amount of mass prescribed by the density model in the 
same angular sector within the radius ri(t) is, 



Am/, = A9 



n(t) 



Ai(ri, 0)ri dn 



oc rj{t)n{n{t),0)A0. 



(3) 

(4) 



To reproduce the model density using regular orbits alone 
we require that 



dlw(I) 



Am/., 



(•») 



where we have taken time averages (as indicated by the 
brackets) of all the times when an orbit comes to the same 
angle, and have summed up all regular families by integrat- 
ing the second integral /. 

It is straightforward to verify this result for the scale- 
free logarithmic disc of K93. In this special case fi oc r~ 1 s(9), 
so Amjj oc ri(t)s(9)A9. Eq. dBf) can then be rewritten in a 
form similar to K93, f dIw(I)(A(t/T) I>: j/ri(t)) oc s(0)A0. 

2.2 Curvature at the major axis 

For the time being we assume that the angular momentum 
of a regular orbit is always non-zero everywhere in the orbit 
(which is clearly not the case near the turning point of a 
boxlet orbit), so that 81 7^ 0. Taking the limit A0 — > 0, we 
have 

^ e ^\\ei{t)\\-\ (6) 

and upon substitution in eq. (g) and (|3j), we obtain 

1 



dlw(I) (r) = const, 



fi(r,0)\\J\\ 



(7) 



where the weights w(I) are non-negative and J = r 2 is the 
angular momentum. 

To examine how the orbits fit the curvature of the den- 
sity near the major axis, we take the double derivative of 
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Figure 1. Various reference orbits (thin solid lines) with the 
same energy in a two-dimensional scale-free elliptic disc with a 
surface density power-law slope 7 = 1 — a = 0.8 and axis ratio 
q = 0.7. The orbits shown are a regular banana orbit, a fish orbit, 
two high resonance orbits, a stochastic orbit and a tube orbit. 
Compare the shape of the orbits with the shapes of the over- 
plotted zero- velocity curve (thin dashed line), the model density 
contour (dash-dotted line) and the contour made by summing up 
the spatial distribution of the 'cousin' orbits, the scaled and re- 
flected copy orbits (the outer solid lines). Near the major axis 
the contours of these orbits bend in the opposite direction of the 
elliptical contours of the model density. The similar contours for 
tubes bend in the right direction, but not strong enough. 



eq. (M) with respect to 8, and substitute the result in the 
equations of motion 



J 



06 



r 3 Or ' " ; 



(8) 



The angular momentum J is nearly constant in the vicinity 
of the major axis where the torque —86/88 — > 0, because the 
force is radial at the symmetry axes. Evaluating the angular 
derivatives at the major axis of the potential (8 — 7r/2), we 
obtain the following simple expression (see Appendix A for 
a derivation): 



Figure 2. As Figure 1, but for an elliptic disc with a surface 
density power-law slope 7 = 1 — a = 1.2 and axis ratio q = 0.7. 
Only the fish orbit supports the potential near the major axis. 



and we have written K r = £ , and Kg = r a , as the radial 

ro r <p ra r q> 

and tangential parts of the kinetic energy scaled by the virial 
rd r 6. The quantities qT and g~ , defined as 



Qx 



1 + 



rdrX 1 



X = [6(r,8),v(r,8)} 



(11) 



describe the curvatures of the equal potential and equal 
density contours at the major axis (alternative expressions 
in terms of double derivatives of s(8) and p(8) are given 
by eqs Af: and A.7). Here qy^,^\ are the axis ratios of the 
best-fitting ellipses to the potential, respectively density, 
contours near the major axis. A similar result can be de- 
rived for any angle 8, but the expression is much simpler 
at the symmetry axes because all first derivatives vanish: 
dg/j,(r,8) — dg6(r,8) = OgJ — 0; the result on the mi- 
nor axis turns out to be not very useful in providing con- 
straints to models. We remark that, e.g., circular orbits 
in an axisymmetric system trivially satisfy eq. (kf) because 
g M = q<f, = A = 1. 



diw(i) { [ q ;- 

where 



(1+7- 



^i) 



X^il + ^Kr + q^Kelg^^ > 0, 



= 0, 



(9) 



(10) 



3 ORBITS IN ELONGATED DISCS 

Consider scale-free discs with surface density of the form: 

(i{r, 8) oc r" 7 (cos" 8 + q n sin" 8)' » . (12) 
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Figure 3. The parameter space for scale-free discs: cusp-slope 7 
versus axis ratio q. No self-consistent elliptic discs can occupy the 
shaded region to the left of the solid line (where 1 — q > 7/2), 
due to the wrong curvature of orbits at the major axis. Regions 
to the left of the dashed curve (h < 0.5, see §5.2) and/or the 
dotted line (C < 0.707, see §5.3) are also restricted by the shapes 
of the available orbits. The crosses indicate logarithmic elliptic 
discs which are found to be likely non-self-consistent in previous 
numerical studies by Kuijken (1993). The error bar indicates the 
cusp strength for the nuclear bar in NGC5121 observed by HST. 



The gravitational potential of these discs is not exactly el- 
liptical, and is formally given by eq. (nh. It can be computed 
with harmonic expansions as in SZ98. When n — 2, the 
surface density is stratified on ellipses of axis ratio q. 

Figures pi and show the density distribution of the 
dominant orbits in two such discs. The heavy solid line in 
each panel indicates the shape of the orbit after all the scaled 
copy orbits and their reflection images are summed up in 
a scale-free fashion; for clarity only the right hand part is 
drawn. We find that near the major axis the density con- 
tributions of all tube orbits, and nearly all boxlet orbits, 
curve in the opposite sense as the model density; the excep- 
tions are the fish orbits in some potentials (Miralda-Escude 
& Schwarzschild 1989; K93), which sometimes support the 
curvature of the model density. 

The above result suggests that a large fraction of scale- 
free elliptic disc potentials cannot be made self-consistent 
near the major axis. While the tube orbits are always 
anti-aligned with the potential (e.g., de Zeeuw, Hunter & 
Schwarzschild 1987), once the 'backbone' of a cored poten- 
tial, namely the set of box orbits, is replaced by the boxlets 
and stochastic orbits of a cusped model, there may be no or- 
bits left to support the potential near the major axis. Since 
a general property of all orbits is that the amplitude of the 
orbital angular momentum always reaches a local maximum 
on the major axis, an orbit tends to spend more time away 
from the major axis than on the axis. As a result, the curva- 
ture of orbits is often opposite to the real density even after 
all 'cousin' orbits are added, which limits the range of non- 
axisymmetric models with cusps that can be self-consistent. 
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Figure 4. The parameter space of non-axisymmetric disc po- 
tentials. Self-consistency is ruled out for a large fraction of discs 
with elliptical density contours (indicated by small crosses) and 
for the whole sequence of cuspy, separable potentials of axis ratio 
0.5 studied by Sridhar & Touma (1997) and Syer & Zhao (1998) 
(figure-eight symbols). The axis ratio here is for the equal poten- 
tial contours, although they are generally not exactly elliptical. 



4 FITTING THE CURVATURE AT THE 
MAJOR AXIS 

We now present a more quantitative analysis of the curva- 
ture of the regular orbits. A simple and necessary condition 
for self-consistency near the major axis of the model follows 
from eq. (Mj: at least some orbits should have a negative 
\ql 2 — (1 + 7 — 7A)] and others a positive value in order to 
place the average at zero. This places a range on the possible 
curvature (or axis ratio) of the potential, as follows: 



[0, 1 +7- 7A max ]max < q$ < [1, 1 + 7 - 7^min]min. (13) 

Since oo > A max > A > A m i n > by definition (cf eq. (pXJ)), 
a clean non-trivial limit is that 



<?<A > 



1 



\/TTT 



(14) 



This condition turns out to be quite powerful. It implies that 
discs with shallow density cusps (small 7) cannot have a 
highly elliptical potential because they would violate the re- 
quirement of the curvature. The curvature condition is triv- 
ially satisfied in the Keplerian regime, since the potential 
becomes spherical {q, 2 — ■» 1). This would be the case near a 
central black hole, or at very large distance of a model with 
a finite mass. 

Our analysis nicely confirms the result of the semi- 
analytical study of SZ98 that all ST97 discs are non-self- 
consistent. Eq. (nTj) is always violated for these discs, be- 
cause 

qj, = - 1 < ■ X for ST97 discs with 7 < 1; (15) 

V3- 7 V 1 + 7 

ST97 models with 1 < 7 < 2 have unphysical negative den- 
sity regions. 

More interesting is the application of the curvature 
criterion (114) to the entire class of elongated discs. Self- 
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consistent elliptic discs (n = 2 in eq. {J12J) ) must have strong 
enough cusps such that 

-2 



1>1* 



2(1 -ft.), 



(16) 



where g M = o is the axis ratio of the surface density. 

Figure H shows that about 50% of the parameter space 
in the cusp strength vs axis-ratio plane of elliptic discs can 
be simply ruled out on the basis of examining the curvature 
at the major axis. Figure shows a similar diagram, but for 
the potential of elliptic discs and the ST97 models. When 
applied to the specific case of the 7 = 1 elliptic disk, the 
lower limit on the allowed axis ratio provided by eq. (114) is 
in harmony with the numerical estimates of K93. 

Perhaps unexpectedly, our criterion shows that non- 
axisymmetric discs with a shallow cusp are easier to rule 
out than those with strong cusps, while the opposite has 
been suggested for three-dimensional models, based on nu- 
merical experiments with two classes of potentials, one with 
a finite force at the center (7 = in projection), the other 
with the force diverge as r" 1 (7 = 1) (Merritt & Fridman 
1996). These two opposing suggestions here may well reflect 
the still incomplete coverage of parameter spaces of both ap- 
proaches: we are restricted to 2D scale-free discs, but we gain 
better coverage of the range of cusp slope 7 and ellipticity 
1 — q owing to our analytical approach. 



5 DISCUSSION AND CONCLUSIONS 

We have shown that many scale-free elliptic discs cannot be 
self-consistent, and we have given specific criteria for pick- 
ing out non-self-consistent models. The major axis is a local 
maximum in terms of the angular momentum of an orbit, 
which often translates to a local minimum in terms of the 
fraction of time spent by regular boxlets and tube orbits. 
This is opposite to that required by the surface density dis- 
tribution. A clean result from this analysis is shown by the 
forbidden zone in Fig. 0, which implies that galactic nuclear 
discs with a shallow density profile (0 < 7 < 0.3) are neces- 
sarily nearly axisymmetric. 



5.1 Hidden assumptions 

Our results apply only to orbits which cross the disc major 
axis with a finite angular momentum ||J|| = r 2 ||#||, so that 
\\0\\~ and its derivatives are well-defined at 9 = 7r/2. This 
condition clearly breaks down for the axial orbits and/or 
the stochastic orbits. The two are perhaps the same since 
the axial orbits are destabilized by the divergent force at 
the center, and become stochastic. 

Stochastic orbits make our method problematic be- 
cause they could reach J = virtually anywhere. And it is 
quite inevitable that they will be populated during galaxy 
formation. However, the way they are populated must be 
restricted in an equilibrium model because stochastic or- 
bits are slowly-evolving orbits. The only sure recipe to mix 
stochastic orbits of the same energy (-Bo) together into a 
time- independent building block (a super-orbit) is with a 
distribution function S(E — Eq) (Zhao 1996); this distri- 
bution mixes in the regular orbits of the same energy Eq 
as a side effect. In our scale-free models these super-orbits 
form one family, with the relative weights of all members 



related by a simple scaling (§2). Fortunately for the present 
problem, the spatial distribution of this super-orbit family 
is always as round as the equal-potential contours. Since the 
equal-density contours are generally more elongated, includ- 
ing the super-orbits does not relax our constraints on curva- 
turcf] So in summary, our result remains valid if stochastic 
and/or axial orbits are populated with a time-independent 
distribution. 

Can a small operation near the major axis of the density 
model bring the model to self-consistency? One possibility 
is to modify the shape of the equal-density contour near 
the major axis by a significant amount (without violating 
the Poisson equation and the positivity of the density) such 
that it is rounder than that of the equal potential contour 
(fti > Q<t>) locally. It is conceivable that by heavily popu- 
lating the super-orbit component, which is stratified on a 
set of mildly elongated equal-potential contours, one might 
counteract the wrong curvature from regular orbits at the 
major axis. Whether self-consistency can be restored this 
way also depends on the maximum fraction of mass which 
can be allocated to the super-orbit component. 



5.2 Expanding the forbidden territory 

Can self-consistent scale-free discs be ruled out on the right- 
hand side of the 'forbidden zone' in Figure fa? Perhaps this 
is likely: if regular orbits always have a finite (rescaled) an- 
gular momentum h when crossing the major axis, and or- 
bits with less angular momentum are stochastic, then the 
minimum angular momentum sets up a barrier to prevent 
any regular orbit from falling to the center or touching the 
zero-velocity curve on the major axis. Here h is defined by 
< h = Jmin/^B < lj where Je is the maximum angular 
momentum allowed for an orbit of energy E. There should 
be a gap in phase space filled by stochastic orbits which 
separates the unstable periodic axial orbit from the regular 
boxlets. Such a stochastic gap is typically seen in the start 
space on the zero- velocity surface of three-dimensional mod- 
els, and in surfaces of section along the major axis of two- 
dimensional models (e.g., fig. 5 of Schwarzschild 1993). The 
minimum angular momentum sets up a dynamical boundary 
to r, rd r 4>, r and A (cf eq. (|lC])) at the upper and lower ends. 
The upper limit on A typically helps little to tighten the con- 
straints, but the finite lower limit of A, namely A m i n (ft) as a 
function of h, propagates to a more stringent upper limit on 
the flattening (cf eq. (113)) and pushes the limit of the forbid- 
den region further to the right in Figure pi The parameter 
space to the left of the dashed line can be ruled out as long as 
the stochastic gap is wide enough (h > 0.5). Typically only 
fish orbits or higher-order resonant orbits can have a small h, 
and we suspect that it is difficult to reach self-consistency 
for much of the region to the right of the forbidden zone 



1 It is possible to modify the shape of a super-orbit by subtracting 
the densities of regular orbits from it — always keeping the den- 
sity everywhere non-negative. It is not clear whether the result- 
ing component (a nearly spherical distribution with many worm- 
holes) will be more elongated than the density model. Even if it 
were to do so, we suspect that it would not relax the constraint 
on the intrinsic shape of the model unless this single component 
is also very heavily populated during galaxy formation. 
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in Fig. H without heavily populating the higher-order res- 
onances. This hypothesis is also supported by the range of 
non-self-consistent models of K93 (indicated by the crosses) , 
which extends to our h = 0.5 line. 



5.3 Constraint from the density contrast 

The density ratio of the minor vs. major axis should also set 
limits on the ellipticity of self-consistent discs. As shown by 
SZ98, if 



C 



Kr,0) 

Mnf) 



0(r,O) 



ts 



(17) 



is defined as a measure for the model density contrast be- 
tween the minor and major axes, where 1 — e M and 1 — e<p 
are the minor-to-major axis ratio of the equal density and 
equal potential contours respectively, then C must have a 
non-trivial lower and upper bound, 



C n 



<c <c„ 



(18) 



set by the shape of orbits in the potential. A model con- 
structed by populating only the thinnest banana orbit in 
a potential cannot be flatter than a certain value because 
of the fact that even the thinnest banana orbits spend a 
fair amount of time near the minor axis. Populating thick 
boxlet orbits, higher-order resonant orbits and tube orbits 
tend to make the shape rounder as these orbits come more 
frequently to the minor axis than the thin banana orbits. 
Figure H also shows the curve with C = 0.707 m l/v2 in 
the q vs. 7 plane. This curve is interesting for reference as it 
corresponds to a configuration where only the thinnest ba- 
nana orbits in a ST97 potential are populated. In contrast 
a C — 1 curve (which means e M = e^) would correspond to 
a configuration where only the f(E) super-orbits are popu- 
lated. The realistic configurations are likely in between these 
extremes. Combined with the criteria from the curvature, 
Figure tfl suggests that about 2/3 of the parameter space of 
cusped elliptic discs are ruled out. 



5.4 Implications for nuclear discs and bars 

Our curvature constraint on the existence of self-consistent 
scale-free elliptic disks can be applied to the properties 
of observed stellar bars and nuclear discs. The dynam- 
ics of these highly flattened systems are essentially two- 
dimensional with the small vertical oscillation completely 
decoupled from the motion in the plane. In the absence of 
a curvature criterion for three-dimensional models, it is pre- 
mature to discuss the parameter space for elliptical galaxies. 
Tumbling stellar bars have scales set by the corotation 
radius. Numerical experiments show that regular boxlet or 
loop orbits have a central 'hole' of finite size compared to 
their apocenter. So one expects a region infinitely close to 
the center where no 'large' orbits will ever visit. The dynam- 
ics of the very nucleus will be dominated by 'small' orbits in 
situ. Such a situation would not be possible for models with 
a finite core. We know of at least one set of strongly elon- 
gated bar models — the rotating Freeman (1966) elliptic 
discs, which can be built self-consistently for any axis ra- 
tio. Our curvature condition clearly does not apply to these 
models with an analytical core (7 = 0) as the finite force at 



the center stabilises axial orbits and box orbits. But it can 
be applied with confidence to the central regions of cusped 
elliptical nuclear bars. And our results show that these can- 
not be strongly elongated; their axis ratio in projected light 
should satisfy the major axis curvature condition eq. wM- 

Our curvature constraint on elongation comes from the 
general property that the orbital angular momentum peaks 
near the major axis, so that the orbital density generally has 
a local minimum on this axis, whereas the model density is 
largest there. We therefore expect it to hold to some extent 
for general discs, as long as a genuine cusp (7 7^ 0) in the 
central surface density creates a divergent force at the center 
and destabilises any box orbit. An outer truncation of the 
mass distribution and the tumbling motion of the bar can 
greatly change the potential at large radius, but at very 
small radii the model reduces to the static self-gravitating 
scale-free case. This happens at radii well-inside the ILR 
of any rotating potential (but well beyond the sphere of 
influence of any small central black hole if it exists in bars) 
such that the centrifugal and Coriolis forces are negligible 
compared to the self-gravity of the cusp. 

Measurements of nuclear cusp slopes for galaxies with 
nuclear bars are as yet scarce. The one object for which the 
light distribution has been measured with HST resolution, 
NGC 5121, shows a steep light profile with 7 = (0.75 ± 
0.05) and a projected axis ratio (0.67 ± 0.05) (Carollo & 
Stiavelli 1998). This data point, as shown in Fig. H, is safely 
outside the forbidden zone with or without correcting for 
the inclination. Imaging of other nuclear bars are needed to 
establish whether there is a zone of avoidance for the shapes 
of cusped galaxies. 

It is a pleasure to acknowledge helpful discussions with Vin- 
cent Icke, Konrad Kuijken, Doug Richstone, Massimo Sti- 
avelli, Scott Tremaine, and Dave Syer, and thoughtful sug- 
gestions by an anonymous referee which helped to improve 
the presentation of our results. 



APPENDIX A: DERIVATION FOR DOUBLE 
DERIVATIVE OF V 

The lengthy calculation of the double derivative of T with 
respect to 6 is greatly simplified near the major axis (9 — §), 
where we can safely ignore any first or third derivatives of an 
even function with respect to time t or angle 6, since these 
tend to zero. 

First we decompose the double derivative of F (cf. eq. M) 
at 9 — ? to that with respect to radius r, to angular mo- 
mentum J, and to the angular part s(0) of the density, 

d 2 F d 2 J dV d 2 s(<9) 

(Al) 



rd6> 2 ' 2 Jd g2 r - ld Q2 sd6 2 ■ 

With the help of the equations of motion (eq. fel) we find 



d 3 J 

Jd6 2 



dV 

r"'d6 2 



= 7 



dJ 
J6d6 

r .2Q2 



" (1+7) bJ 



r(P 



+ 71 



(A2) 
(A3) 
(A4) 
(A5) 
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Now rewrite the second derivatives of the density /j, and the 
potential <j> in terms of the curvatures of the density and 
potential, which by definition (cf. eq. lilt) satisfy 

q ^-i = J|L = «M, (A6) 



rdr/J. s —7 ' 



(A7) 



* rd r (f> p 1 — 7 

Rewriting f 2 and # 2 in terms of K r and -ftTg, we obtain 

jq^lfe* = ^ fc 2 - (1 +7)(1 -7«r) + 7^ 2 ^] ,(A8) 
which immediately reduces to eq. (Mj. 
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